Numerical Analysis of Quasiholes of the Moore-Read Wavefunction 
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We demonstrate numerically that non-Abelian quasihole (qh) excitations of the v = 5/2 fractional 
quantum Hall state have some of the key properties necessary to support quantum computation. 
We find that as the qh spacing is increased, the unitary transformation which describes winding two 
qh's around each other converges exponentially to its asymptotic limit and that the two orthogonal 
wavefunctions describing a system with four qh's become exponentially degenerate. We calculate 
the length scales for these two decays to be £jy « 2.7 £o and £g « 2.3 respectively. Additionally 
we determine which fusion channel is lower in energy when two qh's are brought close together. 



The proposal to use quantum Hall states as a platform 
for qu antum computation has spurred a great deal of 
interesiJS^m. These quantum Hall systems are believed to 
have natural "topological" immunity to decoherence and 
therefore hold particular promise for quantum compu- 
tation. In so-called non-Abelian quantum Hall systems, 
the ground state is highly degenerate in the presence of 
quasiparticles (qp's), and this degenerate space can be 
used to store quantum information. Operations on this 
space are then performed by adiabatically dragging qp's 
around each other, thus "braiding" their world-lines in 
2+1 dimensions. 

Although there is currently no definitive experimental 
evidence that non-Abelian quantum Hall states even ex- 
ist, the community now strongly suspects^ that the quan- 
tum Hall plateau observed at Landau level (LL) filling 
fraction = 5/2 is the non-Abelian Moore-Read (MR) 
phase^ (or its closely related particle-hole conjugate^). 
While the MR phase is, strictly speaking, not capable 
of universal topological quantum computation (compu- 
tation by braiding qp's around each other at large dis- 
tances), a scheme has been devisecPthat in principle al- 
lows error free quantum computation by supplementing 
these topological processes with nontopological processes 
where qp's are moved together and allowed to interact. 
Furthermore, the MR phase is frequently viewed as the 
simplest paradigm of a non-Abelian state of matter, and 
is therefore a logical starting point for detailed analysis^. 

In order for topological (or partially topological) 
schemes for quantum computation to be scalable (i.e., 
to allow large scale quantum computation) , a number of 
crucial conditions must holcP. Condition (1) As all of 
the qp's are moved apart from one another, the split- 
ting of the energy levels of the putatively degenerate 
ground state space must converge to zero at least as fast 
as e~ R ^ B where R is the minimum distance between 
qp's. In the literature, there has been numerical work 
suggesting that condition (1) may not be true^ for the 
MR state. One of the goals of our work is to perform more 
precise numerical calculations to determine whether this 
numerical conclusion holds up to more careful scrutiny. 
Condition (2) As qp's are moved apart from each other, 
the unitary transformation that results from adiabati- 



cally dragging one qp around another must converge to 
its asymptotic limit at least as fast as e~ R ^ u . For the 
MR st ate, s everal theoretical arguments suggest that this 
is tru d 8 * 9 * 10 -* however, in these theoretical works the pre- 
cise length scale £u remains unknown. Presumably £e 
and £jj are both on the scale of a magnetic length mul- 
tiplied by some number of order unity. If this number 
of "order unity" happens to be very large, it could in 
principle start to cause trouble for practical implementa- 
tion of topological schemes. We will explicitly determine 
both £[j and £e numerically. Finally, Condition (3) One 
must be able to measure the topological quantum num- 
ber associated with a group of qp's. Proposals have been 
made that such quantu m numbers can be measured us- 
ing interferometrjEMII] However, this scheme has turned 
out to be very difficult experimentally. Another possi- 
ble way to measure the topological quantum number of, 
say, two qp's, is to move the qp's microscopically close 
and precisely measure the force between them (or equiv- 
alcntly the energy change of moving them) . While this 
may not sound any easier, it nonetheless proposes a dif- 
ferent route to making this measurement should interfer- 
ometry prove to be impossible. In the current paper we 
will attempt to numerically evaluate this energy change 
and show how it reflects the quantum number of a pair 
of qp's. See Ref. [TS] for a similar analysis of the Kitaev 
model. 

Our numerical work is performed on a spherical geome- 
try with a monopole of flux at the center of the sphere 
and N electrons on the surface. For the MR state^, is 
given by = 2N — 3 + n q h /2 and n q h is the number of 
quasiholes (qp's with positive charge). The radius of the 
sphere is (N l p/2) 1 / 2 £o where £q is the magnetic length. 
For the the purpose of stating the decay lengths £j/ and 
£e, the distance R between quasiholes (qh's) will be writ- 
ten in terms of the chord length. The definitions of £e 
and £[/ are given below. It should be noted that while 
their precise values depend on the the particular qh con- 
figurations used in our calculations, alternate definitions 
for different qh configurations will give results that only 
differ by factors of order unity. 

We consider the MR wavefunction (wf) in the pres- 
ence of qh's which is defined as the zero energy space 
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FIG. 1: (color online) Statistical phase for winding one qh 
around another in the spherical geometry as defined in the 
text. Data is shown for the case of even total numbers of 
electrons for which the topological quantum number of the 
pair of qh's is 1. In (a) we have plotted the statistical phase 
versus the chord distance R between the two qh's and fit the 
data using cos(a^ + b) for the oscillatory part and a decay 
term that is either exponential, Gaussian or power law. We 
fit the data starting at R = 6.48 lo and find the value of the 
reduced x 2 ls smallest for exponential decay with a value of 
1.42 while for power law and Gaussian decays it is 7.22 and 
5.52, respectively. The good fit to exponential decay is also 
confirmed when we plot the absolute value of the data on log 
and log-log scales (b) and perform linear fits to the extrema of 
the oscillations. The linear fit is clearly much better on the log 
plot demonstrating that the oscillations decay exponentially 
rather than as a power law. 

of a special short-ranged three-body interaction^. Al- 
though this is just a model interaction, the ground state 
wf turns out to be an accurate approximation for more 
realistic interactions^. Thus our calculations are varia- 
tional in nature. Pairs of qh excitations carry the topo- 
logical quantum number "1" or "-0" which represent the 
two states of a qubit, and the degeneracy^ of a system 
with n q h qh's is 2t x . 

We start by considering the case of two qh's for which 
the ground state is unique. In this case we can address 
condition (2) above by calculating the braiding statistics 
of these two qh's. To do so in the spherical geometry, 
we compute the Berry phase accumulated when one qh 
is moved adiabatically around the equator of the sphere 
while the second qh is held fixed first on the north pole 
and then on the south pole. Both these Berry phases 
have contributions from the statistical phase associated 
with the two qh's, and the Aharanov-Bohm phase due 
to the applied magnetic field. To isolate the statistical 
phase we therefore compute the difference between these 
two phases. In the planar geometry this difference would 
correspond to the change in the Berry phase when one 
qh is moved in a closed loop while a second qh is held 



fixed first inside the loop and then outside the loop. 

The Berry phases are all calculated numerically using 
a Monte Carlo method essentially identical to that de- 
scribed in Ref. HH We use the MR wf with two qh's, 
which is not an exact wf for the realistic Coulomb in- 
teraction, but is quite accurate nonetheless as prior nu- 
merical work has demonstrated^. When we drag qh's 
we can think of having added a highly localized potential 
well to the system whose position moves as a function of 
time. However, since our Berry phase calculation does 
not involve a detailed Hamiltonian per-se, our results are 
independent of the form of this potential well. (Further 
details of the methods used will be given in Ref. [T7)l . For 
the cases of either an even or odd number of electrons 
on the sphere, the two qh's together must have topologi- 
cal quantum numbers 1 or ijj, respectively. The statistical 
phase is then expect ed 1 | 4 | 8 | 1Q | to converge either to zero (if 
the quantum number is 1) or it (if the quantum number is 
ip) as the distance between the qh's is increased. Indeed, 
for an even number of electrons, we show in Fig. [TJthat 
as the sphere is made larger, the convergence is exponen- 
tial and the decay scale is roughly £u ss 2.7 1§. Similar 
results were obtained for the case of an odd number of 
electrons where the phase converges exponentially to 7r 
with roughly the same decay scale. 

The difference between the even and odd case can be 
interpreted as the non-Abelian component - i.e, the part 
of the phase that depends on which topological sector 
the two qp's are in. We conclude that this non-Abelian 
contribution does indeed converge exponentially with in- 
creasing system size as desired by condition (2). (Ide- 
ally we would like to determine the unitary transforma- 
tion that occurs on this two dimensional ground state 
space when particles are braided around each other as 
in Ref. 1161 However, we have found that it is currently 
numerically too demanding to demonstrate exponential 
convergence in this more complicated situation). 

The oscillations in Fig. [I] (and in the later Figures) are 
not unexpected. In the closely related system of a p- 
wave paired superfluid, the oscillatin g form of the wave- 
functions can be calculated explicitl y 18 * 19 !. However, in 
this quantum Hall system those results would only be 
qualitative. 

To address condition (1) above, we now turn to the 
case of four qh's and restrict ourselves to an even num- 
ber of electrons. We implement a trial wf approach using 
the Moore-Read wf with qh's, which is the ground state 
of a special three body interaction, but we will evaluate 
its energy with a 1st excited LL Coulomb interactional, 
for which the MR wf is not the exact ground state. We 
nonetheless expect this hybrid approach to give accurate 
results because the MR wf is an extremely accurate ap- 
proximation of the exact ground state of the Coulomb 
interaction. Our calculation is also an exact statement 
about the lowest order perturbation of the special three 
body interaction towards the Coulomb pointful 

For the MR wf with four qh's, there are two puta- 
tively degenerate ground state wavefunctions^H Using 
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a spherical geometry, we place four qh's on the corners 
of an equilateral tetrahedron and implement standard 
Monte Carlo procedures to evaluate the energy splitting 
between the two eigenstates of the interaction within the 
two dimensional ground state subspace (Details will be 
presented in Ref. [17j). Results are presented in Fig. [2] 
as a function of system size, and indeed it appears that 
the two blocks become degenerate exponentially as the 
distance between qp's increases as required by condition 
(2) with a decay length of £e m 2.3 £q. This result ap- 
pears to contradict results of Ref. |7] which claimed an 
algebraic rather than exponential decay. Although the 
methods used by the two works are essentially identical, 
in Ref. 7] the MR state is studied in the lowest LL (where 
MR is known not to be a good trial wf) whereas we have 
studied it in the first LL where it is known to be a very 
good trial wf and thought to be experimentally relevant. 
Differences with Ref. |7] could occur because of differing 
levels of Monte Carlo error as well. The numerical diffi- 
culty of collecting data is substantial, so admittedly our 
error bars are currently somewhat larger than desirable. 
However, we will continue to collect data and these re- 
sults will almost certainly improve. 
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FIG. 2: (color online) Energy splitting of the eigenstates on 
a sphere with four qh's as a function of system size. The four 
qh's are placed on the corners of an equilateral tetrahedron. 
The data is fit with the function cos(ay/N + b) multiplied by 
an exponential, Gaussian, or power law decay function. The 
distance between particles grows as yN, so this is essentially 
the same fit as used in Fig. [I] The data is fit starting with 
N = 12 (not shown because the fits are nearly identical at the 
small N values). The reduced \ 2 values for the exponential, 
Gaussian, and power law fits are 3.39, 7.73, and 6.49 respec- 
tively, which helps confirm our expectation that two energies 
become exponentially degenerate as the qh's move apart. We 
have not shown the log and log-log plots here because such 
plots discard the sign, and in the absence of a higher density 
of points, are hard to interpret. 

Finally we turn to the issue of measurement, condi- 
tion (3). Here, we start with four qh's at the corners 
of a tetrahedron on a relatively large sphere (N = 40) 
where the two ground state wavefunctions are close to 
degenerate. We then move qh 2 close to 1 and observe 
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FIG. 3: (color online) This figure shows the total energy of 
the N = 40 system with four qh's as qh 2 is moved closer to 
its pair, qh 1. Qh 1 is located at the north pole and qh 2 is 
moved along a path that keeps a certain analytic forrrP for 
the trial states |1) and \ip) precisely orthogonal. The two 
states of the system, |1) and which are nearly degenerate 
when the qh's are well separated, split as the qh's approach 
each other. The inset shows the energy needed to move qh 
2 within Iq of qh 1 for different values of TV. We find that 
it takes more energy to bring the two qh's together when the 
system is in state |1) than in state and that the energy 
splitting is on the order of 0.01e 2 /e£ . 



the change in energy of the two wavefunctions. It turns 
out (and we will show in detail in Ref. T?j that if we 
choose to move the qh's together along an appropriately 
chosen path, then the conformal block wavefunctions de- 
fined in Ref. [5| diagonalize the interaction. These two 
conformal block wavefunctions, known as |1) and |^>), 
are constructed such that the pair of qh's 1 and 2 have 
topological quantum number 1 and ip respectively. The 
results of such a calculation are shown in Fig. [3] We see 
that the energy of moving two qh's together is always 
positive simply due to the Coulomb repulsion. However, 
the energy is substantially greater when the two qh's are 
in the |1) state compared to the \ifj) state. To our knowl- 
edge, this result was not predicted and may be attributed 
to the fact that the electron density vanishes in the |1) 
state when two qh's approach each other, but remains 
non-zero in the \ip) state resulting in a more extended 
objectP. Thus our calculation makes the first mapping 
between a proposed measurement of the energy of two 
qh's and what this would indicate in terms of determin- 
ing their topological quantum number. (Obviously if one 
were moving a qp together with a qh, the |1) state would 
have lower energy). We also point out that this result 
may be significant in regards to the possibility of qh's 
condensing into daughter states as proposed in Refs. |2"T1 
and El 

The magnitude of the energy splitting of the two states, 
when two qh's are very close to each other (within £q), 
is measured to be roughly O.Ole 2 /c£q which in a real sys- 
tem corresponds to roughly 1 K, a rather small energy 
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to be measured. To make matters worse, this measured 
energy should be considered to be an upper bound, as 
mixing with states above the gap will be substantial and 
could easily reduce this energy scale (the experimentally 
measured gap itself is less than 1 K in the very best 
samples, although theoretically without disorder the gap 
could be almost 2.5 K. See Ref. |T|and therein). Nonethe- 
less, this numerical work gives the first order of magni- 
tude estimate for how large the splitting due to topolog- 
ical quantum numbers is likely to be compared to the 
overall Coulomb energy between the two qh's. 

The decay length scales and energy scales that we cal- 
culate above are also extremely relevant to majorana tun- 
neling. Plugging in real numbers, we find that at a sep- 
aration of about 0.1 micron, the energy diff erence will 
be about 80mK with the sign of the tunneling ampli- 
tude depending sensitively on the distance. These scales 
have implication to interferometry experiments^ where 
tunneling occurs between edge and bulk as well as for ma- 
jorana hopping problem where tunneling occurs between 
many bulk qh'a^l 

To summarize, we have used Monte Carlo techniques to 
examine several key properties of the MR wf with qh's. 
Note that because our calculations do not incorporate 



LL mixing terms (which are expected to be small), they 
are equally applicable to the recently proposed AntiP- 
faffian wP. We find that both the unitary transforma- 
tion associated with adiabatic transport and the energy 
splitting of putativcly degenerate states converge expo- 
nentially with increasing distance between qp's, and we 
explicitly extract the decay lengths. Encouragingly, the 
decay lengths are on the order of a magnetic length which 
suggests that qp spacing should not be a barrier to physi- 
cal implementation of topological operations. Further we 
examine the energy splitting that occurs when two qh's 
are moved together. We find that the |1) state of these 
two particles is of higher energy and we measure this en- 
ergy splitting between |1) and Although this energy 
splitting is small, it gives experimentalists another way 
to measure topological quantum numbers in these sys- 
tems. Many more details of this work will be presented 
in an upcoming publication^. 
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